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,-H Abstract 

I"""! The first generation of ground-based interferometric gravitational wave de- 

Q^ tectors, LIGO, GEO and Virgo, liave operated and taken data at tlieir design 

^ sensitivities over tlie last few years. The data has been examined for the pres- 

, "^ ence of gravitational wave signals. Presented here is a comprehensive review 

of the most significant results. The network of detectors is currently being 

^^ upgraded and extended, providing a large likelihood for observations. These 

"Tj- future prospects will also be discussed. 
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1. Introduction 



o 
en 

• ^ GW observations promise to bring new insight in many phenomena, prob- 

rN ing compact objects such as neutron stars and black holes, cosmological mod- 

c3 els of the early Universe and testing general relativity in the strong field 

regime [H 121 El S] ■ With the construction and the successful operation of the 
first generation of long arm laser interferometers in the United States (LIGO), 
Italy (Virgo) and Germany (GEO) [3, EJ [7j the last two decades have been 
incredibly important for gravitational wave (GW) searches in the 50 — 1000 
Hz range. The upgraded LIGO, GEO and Virgo detectors [8l|9], expected to 
come online after 2015, along with the Large-scale Cryogenic Gravitational 
wave Telescope (KAGRA) project in Japan [10] will probe a volume of space 
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1000 times larger than the first generation of detectors yielding a high proba- 
bility for multiple observations. For instance, 40 binary neutron star inspirals 
per year are expected to be observed with the advanced ground-based detec- 
tors [11] at design sensitivity. However, the first GW detection might be in 
the 10~^ — 10"*^ Hz range and come from the analysis of the timing residuals 
for sets of radio pulsars [12]. In the next decade we also hope to see a GW 
space mission fly, accessing the incredible number of sources in the 10~^ Hz 
- 0.1 Hz bandwidth. 

GW observations will complement the electromagnetic (EM) ones: While 
EM radiation is scattered, absorbed and delayed, due to their weak inter- 
action, GWs carry information from the entirety (especially including the 
onset) of cataclysmic astrophysical events. Moreover, because of their emis- 
sion mechanism, GWs carry information about the bulk motion of matter 
rather than on the differential distribution of charges of the source. Asso- 
ciating Gamma Ray Bursts (GRB) to the gravitational signature of binary 
neutron star and black hole mergers or supernova core collapses will solve 
an enigma of more than 40 years [131 El [IS]- After many GW detections 
statistical population studies of the masses of compact binary inspiral sys- 
tems will address many questions pertaining what is now just an hypothetical 
population of stellar mass (up to ~ 10 Mq) and intermediate mass (up to ~ 
1000 Mq) black holes. We also expect to observe a number of GW sources 
that are EM-invisible: for instance for rapidly spinning neutron stars there 
exist various channels through which they could to be emitting GWs while 
appearing EM silent to us. The detection of such radiation could bring a 
wealth of information about the underlying emission mechanism involving 
crustal deformations, magnetic fields, accretion, and the equation of state of 
the neutron star interior [TB] . 

Individually, each ground-based GW detector is unable to pin down the 
location in the sky of the source of a short GW signal. A set of widely 
separated interferometers can determine the source position using timing, 
amplitude and polarization information of the signal seen coherently in the 
network [HI [IB]. For this reason since 2007 the LSC (LIGO and GEO) and 
Virgo have been pooling their data and carrying out their searches jointly. 
In 2007 the LSC was about to conclude its 5^^ science (S5) run while Virgo 
had just started its first science run (VSRl). Since then, another one year 
long joint science run took place in 2009-2010 (S6, VSR2 and VSR3). In late 
2010 the LIGO detectors shut down for advanced LIGO installation. Virgo 
continued commissioning through late 2011 and improved its low frequency 



sensitivity (< 60 Hz) by a factor of 3. A joint run with the GEO detector 
was conducted in summer of 2011 (S6 for GEO and VSR4 run for Virgo). In 
the following sections, we present the most important results obtained with 
the LIGO, GEO and Virgo data taken during these last years in pursuit of 
the main sources of GWs. We will conclude with a brief overview of the 
astrophysical potential of the network of second generation detectors. 

2. Transient sources and multi-messenger analysis 

Short-lived, or transient, GW signals are thought to be generated in a 
variety of circumstances. In some cases, for instance the gravitational collapse 
of a massive star, or the coalescence of two neutron stars, GW emission 
is associated with EM radiation or neutrino radiation. In other cases the 
scenarios are more speculative and refer to objects that have never been 
observed as, for example, in the case of stellar mass black hole in binary 
system or cosmic strings. But all these scenarios could generate transient 
GWs lasting from a few milliseconds up to a few minutes or days in the most 
extreme cases. 

2.1. Inspiral phase of compact binary systems 

One of the most promising GW sources for interferometric gravitational 
detectors is the coalescence of two compact objects: binary neutron stars 
(BNSs) , binary black holes (BBHs), or a binary neutron star-black hole (NS- 
BH) system. As observed with the binary system PSR 1913-1-16 discovered 
in 1974 [19], GWs are emitted as the two objects slowly spiral towards each 
other. They will eventually merge and form a black hole. 

While binary systems spiral for billions of years, only the signal from 
their last instants is detectable: a few minutes to a few seconds before the 
merger, when the GW frequency is sufficiently high to fall in the detectors' 
observational bandwidth. The inspiral phase is well described by analytical 
post-Newtonian theory until the system reaches the innermost stable circular 
orbit. Beyond this orbit strong gravitational forces dominate and the two 
bodies merge in a fraction of a second, forming a single black hole that reaches 
equilibrium through quasi-normal mode GW emission (ring-down phase). A 
BNS inspiral signal will have many cycles in the frequency bandwidth of the 
LIGO and Virgo detectors. For systems with total mass > 100 Mq the inner 
most stable circular orbit is lower than the lowest frequency of LIGO and 
Virgo detectors and only the merger and ring-down phases are visible by these 



detectors for systems of total mass up to 450 M©. Note, however, that the 
existence of intermediate mass black holes (systems with total mass between 
100 and 450 M0) is still under debate. Globular clusters are suspected to be 
the best environment for direct collapse (stalled supernova) of very massive 
first generation stars (population III stars) [20] or for runaway merging of 
young stars pi] leading to intermediate mass black holes formation. 
The GW inspiral signal depends on the physical parameters of the binary 
system: masses, spins and eccentricity. For isolated systems it is generally 
admitted that the binary systems' orbit has had time to become circular be- 
fore entering the detector's sensitivity band, so eccentricity is neglected. Spin 
is also neglected for low mass binaries, which we expect not to affect detection 
efficiency. Under these assumptions rather accurate waveforms can be con- 
structed for signals from the inspiral phase of compact binaries p2] and used 
in matched filtering based searches. Only recently the merger phase of some 
binary configurations has been numerically simulated (see [23| [2^ 125] and 
references therein). A number of waveforms including inspiral, merger and 
ring-down phases, without spin have been constructed and used for searching 
for binary systems with total mass up to 100 and 450 MqI^H] [27]. However, 
these waveforms are subject to large uncertainties for mass ratios > 1:4. Ne- 
glecting spin also introduces losses in detection efficiency for intermediate 
mass black holes searches, and despite the tremendous efforts of the numer- 
ical relativity community, template-banks of fully spinning waveforms with 
spin precession for any mass ratio are not yet available. All in all, matched fil- 
tering searches are limited to systems of total mass smaller than 100 Mq and 
the higher mass parameter space is covered by un-modelled generic transient 
searches focusing on the merger and ring-down phases (~ 1 s long signals). 
The identification of the source parameters (sky location, masses, spins, 
polarization, orientation, ...) will be the step immediately following a GW 
detection. Markov Chain Monte Carlo or Nested Sampling techniques have 
been exploited to do this [281 129] exploiting waveforms with fully precessing 
spins and the most recent theoretical investigations 



2.2. Sources of hursts of GWs 

As seen in the previous section, as the total mass of the compact binary 
system increases, the inspiral part of the GW signal happens at lower and 
lower frequencies leaving only the merger and ring-down phase in the sen- 
sitive band of the instruments. The merger phase is rapid (seconds or less) 
and despite the advances in numerical relativity it is not well modeled. For 



this reason intermediate mass binary black hole mergers are also searched 
using methods that do not rely on a specific waveform model, but rather 
on excess-power techniques whose detection efficiencies do not depend on 
specific matching of wave-shapes. This type of methods were originally de- 
veloped in the broader context of searching for generic short-lived signals. 
The spectrum of transient searches includes searches associated with EM or 
neutrino counterparts as well as broad all-sky, all-times searches: Unforeseen 
transient phenomena involving, say, newly formed black holes or neutron 
stars accreting matter in highly magnetized media could well be associated 
with the emission of bursts of gravitational radiation. Known mechanisms 
for the emission of GW bursts include the gravitational core collapse of a 
massive star, neutron star glitches and cosmic strings. In the remainder of 
this section we will briefiy review these mechanisms. 

The different mechanisms involved in stellar collapse have been studied 
for more than 40 years, and GW emission is predicted in different phases: it is 
anticipated during the collapse itself through the aspherical bulk mass motion 
(core bounce phase). In the post-bounce phase, convection instabilities inside 
and above the proto-neutron star can develop, driven by entropy gradients, 
neutrino heating and standing accretion shock instability (SASI) mechanisms 
(see [311 E2] and references therein). These phenomena may produce GWs 
lasting up to few seconds. After a supernova explosion or a collapsar (black 
hole formation) a certain fraction of the ejected material can fall back, heating 
the proto-neutron star and exciting oscillation modes that generate GWs (g- 
mode pulsations and r-mode instabilities). Furthermore, when a black hole is 
formed, accretion-induced ringing may generate GWs. Current state-of-the- 
art numerical simulations take into account many non axisymmetric effects in 
3D, attempt to model micro-physics (neutrino transport, electron capture, 
realistic equation of state) and include rotational dynamics and magnetic 
fields. Still, we are quite far from having a complete catalogue of all possible 
waveforms. For instance, GWs emitted during the core bounce exhibit a 
sharp peak in the waveform for most of the initial conditions and energy up 
to a few 10~* MqC^ can be emitted [33], but the exact signal shape depends 
on the rotational profile and the nuclear equation of state. Consequently, for 
the searches, we consider only the generic broad features of the waveforms 
in the literature. 

Some pulsars occasionally show sudden discontinuities (glitches) in their 
EM activity followed by gradual recoveries that may last days or years. Soft 
Gamma-ray Repeater (SGR) fiares [3l], anomalous X-ray pulsars (AXPs) as 



well as radio glitches are examples of this and are thought to be consequences 
of star-quakes or of angular momentum transfer between a solid crust and 
a more rapidly rotating loose component of the neutron star interior |35] . 
During such glitches normal modes (mainly f-modes) could be excited and 
generate GWs. 

Cosmic strings predicted in symmetry-breaking phase transitions [3S] or 
cosmic super-strings in string theory inflation models form a dynamical net- 
work. When cosmic strings interact, they generate loops whose oscillations 
can lead to cusps (sharp bends) which move at the speed of light and gene- 
rate burst of GWs [37]. The GW amplitude depends on the tension of the 
string. The (super-) string density depends on the the loop size, the string 
tension and an additional parameter, the reconnection probability, which in 
the case of super-strings is smaller than 1. 

2.3. Multi-messenger analyses 

Knowing the location of a putative source and the time of occurrence of 
a GW event thanks to the existence of an external, EM or neutrino, trigger 
restricts the GW search parameter space and improves the sensitivity of trig- 
gered searches by roughly a factor 2 with respect to all-sky, all-time searches. 
GRBs, long and short, magnetar flares (SGRs and AXPs), pulsar glitches and 
high energy neutrinos events have been used as triggers for short-lived-GW 
searches. In the remainder of this section we will briefly describe the physical 
process underling each external trigger and describe how it could be associ- 
ated with GW production. We will conclude by considering multi-messenger 
analyses in which the GW search triggers follow-up observations rather than 
being triggered by an EM or neutrino event. 

GRBs are the brightest EM phenomena in the Universe and are likely 
associated to cataclysmic stellar events with a typical energy release of 10^^ 
ergs. Long GRBs are related to the core collapse of massive stars [38] 
|39] which, we have seen in the previous section, can generate GWs. Short 
GRBs are thought to be generated when a BNS or a NS-BH merge forming 
a black hole and a dense accretion disk jlHl E]. As we have seen in the 
previous section BNS and NS-BH mergers are among the prime targets for 
LIGO/Virgo. A small fraction of short GRBs may come from SGR giant 
flares |12]. GWs emitted by galactic SGRs could be detected by the first 
generation of interferometric detectors. 

Single high energy neutrino (HEN) events, although not yet directly ob- 
served, could well be emitted in processes that also produce GWs. Since the 
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HEN and GW detectors have uncorrelated backgrounds, requiring a coinci- 
dence in time and source location between the HEN and GW events could 
significantly reduce the false alarm rate of a joint observation and hence 
improve the sensitivity of the search |13]. Plausible scenarios for joint HEN- 
GW emission include again GRBs, as well as micro-quasars and magnetars 



In core collapse supernova, EM waves are emitted in optically thin regions 
far from the core, which provide limited information about the early instants 
of the core collapse. Low energy neutrinos and GWs are emitted deep inside 
the core and would provide insight in supernova mechanism. Moreover, the 
coincident detection of GWs and low energy neutrinos in the absence of any 
EM counterpart would provide evidence for the existence of so called failed 
supernovae, which could occur at a rate comparable with that of visible 
supernovae |36] . 

Unlike EM telescopes that can only observe some part of the sky at any 
given time, GW detectors see the whole sky all the timq^ This means that 
GW observations could be used to alert EM telescopes of interesting loca- 
tions to point: The putative source location of GW triggers identified in 
close to real time could be passed on to telescopes for follow-up observa- 
tions. Being able to measure the full light curve of a supernova or GRB 
prompt emission and afterglow in all wavelengths would help understand all 
the mechanisms involved in the cataclysmic event, including an investigation 
of the host galaxy at the time of the event. Since directional estimates are 
still challenging for a 3 GW detector network (error box of up to tens of 
square degrees for signals just above detection threshold [SIIIH]), large field 
of view telescopes will be required for this during the first GW observations. 
In 2009-2010, a first campaign of ~ 3 months with an heterogeneous set 
of partners (optical telescopes. X-ray and radio instruments) brought much 
valuable insight for the preparation of advanced detectors era GW alert sys- 
tem [l9l[50], pointing out especially the loose accuracy of the GW source sky 
localization. 

2.4- Transient search results 

To date, no GW transient signal has been discovered in LIGO, GEO or 
Virgo data. As a consequence, upper limits on the rate of GW events or on 



^this is strictly true only with a network of GW detectors spread over the world, where 
blind sky regions of one detector can be covered by other detectors 



the amount of energy released have been derived which can be compared to 
astrophysical predictions. 

Predictions of the rate of coalescences of compact binary systems are 
rather uncertain, spanning 1 or 2 orders of magnitude in each direction from 
the most likely prediction. They are based on observations of binary neu- 
tron stars in our galaxy or on simulations (population synthesis) of systems 
involving black holes, which have not been observed yet. The rates are 
then extrapolated to our local Universe considering the stellar birth rate in 
nearby spiral galaxies. Following |llj we will use the following figures as re- 
alistic predictions for the number of expected coalescences per Mpc^ and per 
year for various binary systems (also shown in Figure Q: 10~^ Mpc^^ yr~^ 
for BNS systems, 3 x 10"^ Mpc"^ yr~^ for a 10 M© NS-BH system and 
5 X 10-9 Mpc-=^ yr"^ for a 10 - 10 Mq BBH system. 

The most recent data (from the S6,VSR2-3 runs) have been searched for 
signals from compact binary coalescences of systems with total mass up to 25 
M0. The waveforms used did not include spin and were terminated just be- 
fore the two objects merge. The typical distance up to which the search could 
have detected these systems was 40 Mpc, 80 Mpc and 90 Mpc for the BNS, 
NS-BH and BBH systems respectively. Because no detection was made upper 
limits on the rates of such astrophysical events were placed at the 90% confi- 
dence level: 1.3 x 10"^ Mpc^^ yr"^ for BNS systems, 3.1 x 10"^ Mpc^^ yr-i 
for a NS-BH system and 6.4 x 10"^ Mpc"^ yr"^ for a BBH system as shown 
in Figure [I] [51j. These numbers fold in as priors the results from previous 
runs. The binary neutron star rate upper limit is still two orders of mag- 
nitude higher than what is likely expected. The detection efficiency of the 
search is lower for spinning systems resulting in upper limits on the rate of 
NS-BHs and BBHs that are ~ 16% larger than for non spinning systems. 
Spin is not expected to significantly change the waveform of BNS systems 
visible by the instruments. 

For binary systems with a total mass higher than 25 Mq (BBH system 
regime) the merger phase contributes a non negligible fraction of the total 
signal-to-noise ratio. For this reason in the search [25] for systems with to- 
tal mass between 25 M© and 100 M© the template-bank utilizes waveforms 
that cover both the inspiral, merger and ring-down phases of the signal. The 
resulting maximum distance reach on S6,VSR2/VSR3 data is 300 Mpc for a 
20 - 20 Mq system, which is significantly higher than for any system in the 
low-mass search reported above [HI]. In the search post-Newtonian inspiral 
waveforms are connected to the merger and ring-down waveforms with an ef- 
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Figure 1: Comparison of upper limit rates for BNS, NS-BH and BBH systems [51 . The 
dark gray regions show the most recent upper Umits obtained with LIGO and Virgo data. 
The new hmits are up to a factor of 1.4 improvement over the previous results (light grey 
regions). The lower (blue) regions show the spread in the astrophysically predicted rates, 
with the dashed-black lines showing the "realistic" estimates [TT] . 



fective one-body approach tuned on numerical relativity waveforms [521 153] . 
Spin effects are not included in the template waveforms. Including spin in 
the search in the most general way is well outside what is currently possible, 
due mainly to the lack of numerical relativity simulations over the broad 
parameter space of possible configurations. There exist, however, waveform 
models for spin configurations aligned or anti-aligned with the orbital an- 
gular momentum, at least for mass ratios between 1 and 4. Figure 7 of 
[26j shows that the average reach of the search for such systems varies by 
about 10% with respect to the non-spinning population. Upper limits on 
the rate of coalescences of such systems were derived as a function of the 
component masses. For example, for systems of SOMq-I-SOMq the rate up- 
per limit is 0.7 X 10^^ Mpc~^ yr~^ and the distance reach well exceeds 200 
Mpc. We refer the reader to Figure 5 of ^26j for the complete set of re- 
sults over the different mass bins. Expectations for the rate of coalescences 
of high mass systems are rather uncertain. However there exist predictions 



n for the rate of 3.6l2'g x 10~^ Mpc^^ yr~^ for stellar mass binary systems 
[5^ [55] . In the corresponding mass bin the GW search constraints the rate 
to 5.9 X 10~^ Mpc~^ yr~^ which already includes part of the predicted rate- 
space, but spin effects were not considered in the search. 

For binary systems of higher mass the only part of the coalescence wave- 
form in the instruments' band is the merger and ring-down. A simple ring- 
down search carried out in 2004 LIGO data yielded an event rate upper limit 
of 3.2 X 10~^ Mpc~^ yr~^ for binary total mass between 100 and 400 M0. In 
S5,VSR1 data a search for intermediate mass black hole systems in the total 
mass range 100 and 450 Mq and with mass ratio up to 1:4 was performed. 
This search is an excess power search, it does not rely on a specific wave- 
form model and coherently combines the data from the different detectors 
|27j . The largest effective range, 241 Mpc, is achieved for a 88 -|- 88 M© 
system and in this mass bin the rate upper limit is 1.3 x 10~^ Mpc~^ yr~^. 
The upper limit on the event rate, averaged over all considered masses, is 
9 X 10~^ Mpc~'^ yr^^ which is just an order of magnitude higher than the 
most constraining rate upper limit. As globular clusters (GC) are the most 
likely hosts of intermediate mass black holes, the event rate can be converted 
into a limit on the astrophysical source density: 3 x 10~^ GC~^ yr~^. This 
is still 3 order of magnitudes higher than the astrophysical predictions [11] , 
which is anyway rather uncertain as intermediate mass black holes formation 
is still under debate. 

The most sensitive all-sky generic transient signal search to date was 
carried out on the S6,VSR2-3 data and would have been able to detect with 
50% efficiency a burst (<ls long) of 2.2 x 10~^ MqC^ of gravitational radiation 
with significant energy content around 150 Hz at the distance of 10 kpc jHE] ■ 
The total observational time for this search amounts to 1.74 yr. Upper limits 
on the rate of different types (waveforms) of events in the frequency range 
between 64 Hz to 5 kHz, as a function of signal amplitude at the Earth are 
derived. For signals strong enough (root-sum-square strain amplitude h^ss ~ 
several 10^^° ^^/rr ) the rate upper limits all tend to 1.4 yr^^. The results 
can also be interpreted as limits on the rate density of GW bursts assuming 
a standard candle source isotropically distributee^ For instance for sources 



^based on two observed tight binaries with individual masses of ~ 15 M0 and ^ 20 
•^this is a rough estimate when the distances are within our local universe. 
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emitting 10~^ Mqc"^ around 150 Hz, the rate is bound to 6x 10^'^ Mpc^^ yr^^- 
As explained above, short GRBs are probably associated with the coales- 
cence of NS-NS or NS-BH systems. On the other hand long GRBs have been 
associated with massive star core collapse. Even though only for ~ 10% of 
the observed GRBs there exist redshift measurements it is believed that the 
typical distance for GRB progenitors be on the Gpc scale, well beyond the 
reach of the first generation of GW detectors. However based on observa- 
tions of a few GRBs at much closer distances (e.g. GRB 980425 at 36 Mpc) 
and on indications on the existence of a class of sub-luminous long GRBs, 
one cannot exclude the chance that there might be a detectable gravitational 
signal in association with a GRB trigger. For this reason searches for GWs 
from binary inspirals and searches for generic GW transients have been car- 
ried out specifically around the time of occurrence of 350 GRBs observed 
by gamma-ray satelhte detectors since 2003 |571 |5H1 |59l |60] . The Swift and 
Fermi satellites have been the primary sources of the most recent GRB trig- 
gers targeted with LIGO- Virgo dataQ taken when at least two GW detectors 
were operational. No GW detection resulted from any of these searches. The 
most recent search [57] excludes the existence of progenitors within a certain 
distance from Earth, which is different for every GRB, assuming an emission 
of 10~^ MqC^ at 150 Hz in GWs. The median distance for this set of GRBs 
is 17 Mpc. From the null results of the binary coalescence signal searches the 
corresponding median distance values are 16 Mpc and 28 Mpc for NS-NS and 
NS-BH systems respectively. Further, if one assumes all GRB progenitors to 
emit the same amount of GWs, the set of GRBs can also be considered as 
a whole. Constraints on the cumulative number of GRBs as function of the 
redshift are shown in Figure |2] where a constant co-moving density rate of 
GRBs up to a cosmological distance beyond which all GRBs are at infinity 
is assumed. Extrapolations at the foreseen advanced GW detectors sensitiv- 
ities reach redshift values in the range of current short GRBs observations. 
For long GRBs, one needs to assume a rather optimistic GW emission (10~^ 
MqC^) in order to reach observed redshifts. 

GW observations were able to exclude certain galaxy associations for two 
GRBs, GRB 051103 [61] and GRB 070201 [62], at least if these GRBs were 
due to BNS mergers. The position obtained by the interplanetary network 



"^the list of the analyzed GRBs, their origins and their main characteristics are given 
in publications [57l [58 l [59 l l60] . 
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Figure 2: Cumulative distribution 90% confidence level exclusion as function of redshift 
from the analysis of 151 long GRBs (left) and 26 short GRBs (right) for different hypothesis 
on the GRB progenitors (smooth curves) . The dashed curves are extrapolations assuming 
a factor 10 in sensitivity is brought by advanced LIGO- Virgo detectors. The red stairs 
show the cumulative distribution of measured redhsifts for GRBs [5 7) . 



of satellites for GRB 051103 had a 3 — a error region overlapping with the 
outer disk of the M81 galaxy (3.63 Mpc). GRB 070201's position was 1.1 
degree away from the center of Andromeda-M31 (770 kpc). A GW binary 
coalescence signal from either M81 or M31 would have been well within the 
sensitivity reach of the operating LIGO detectors. The absence of a GW 
detection following searches for such signals allowed to exclude with high 
confidence the association between the GRBs and a BNS merger in their 
apparent host galaxies. However it was not possible to exclude SGRs in the 
apparent host galaxies as the progenitors of these observed short GRBs (~ 
15% of the short GRBs are accounted for as SGRs [12]). 

GW signals around the times of observed X-rays emission of 6 galactic 
SGRs and anomalous X-ray pulsars have also been searched. The amount of 
GW energy to be expected and the mechanism of emission are still poorly 
understood. However, some theoretical models predict a GW energy emis- 
sion as large as 10^^ — 10^^ ergs [US]- Among all magnetars that have emitted 
x-rays when LIGO, Virgo and GEO were taking data, SGR 1806—20 was the 
most energetic transient observed in our galaxy (~ 4 x 10^^ ergs). Given the 
astrophysical uncertainties, no specific waveforms were assumed in the search, 
and 90% upper limits on the amount of GW energy emitted by these objects 
have been set considering white noise bursts with frequencies ~ 100 Hz as 
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representative of GW signals from low frequency oscillation modes and con- 
sidering damped sinusoids with frequencies > 1 kHz as representative of GW 
signals from f-modes. The most constraining upper limit has been obtained 
for SGR 0501—4516 (1 kpc): it ranges from 3 x 10"^^ ergs to 1 x 10^^ ergs for 
un-modelled low frequency GW emission and for f-mode damped sinusoid 
waveforms with 1090 Hz central frequency, respectively [MIES]- It is inter- 
esting to note that values around 10^^ ergs are comparable to the EM energy 
of giant fiareqj Furthermore the f-mode limit of 1 x 10^^ ergs is in the range 
predicted by [63] . 

The Vela pulsar is one of the pulsars with the least stable spin phase 
evolution: Superimposed on the steady decrease in the rotation velocity of 
this neutron star, astronomers observe a few times a year sudden increases in 
its spin. This could be a consequence of star-quakes that adjust the crust of 
the neutron star or rearrangement of the interior material, and are expected 
to excite normal mode oscillations that may emit GWs. On August 12*^ 
2006 a glitch amounting to a relative increase of the rotational frequency of 
2.62 X 10~^ was measured. This glitch happened when the two collocated 
LIGO detectors at Hanford were taking data. A search was performed for 
ring-down GW waveforms from the quadrupolar (spherical harmonic order 
i = 2) f-mode oscillations with frequencies between 1 and 3 kHz and damping 
times between 50 and 500 ms. No GW signal was observed and 90% con- 
fidence level upper limit on the peak strain amplitude of ring-down signals 
have been set between 6.3 x 10~^^ and 1.4 x 10~^°, depending on which of 
the 2i + l harmonics is excited. This peak strain value corresponds to a GW 
energy range of 5 x 10^^ ergs to 1.3 x 10^^ ergs |66j . 

High energy neutrinos (HENs) may be emitted in various circumstances 
in association with GWs and although many scenarios, such as GRBs and 
SGRs, also involve EM emission that is already used as trigger for targeted 
GW searches, it is not excluded that there might be joint GW-HEN emission 
with scarcely or not at all visible EM counterpart (for example associated 
with chocked GRBs). The first GW search triggered by HEN candidates was 
carried out using HEN triggers from the Antares experiment on S5/VSR1 
GW data [67] . The astrophysical significance of the null result that ensued is 
limited, however this is the first ever GW-HEN search and it paves the way 
for searches in the advanced GW detector era when, with GW sensitivity 



''it is however unknown whether EM and GW energies are correlated. 
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increases of up to a factor of 10 it is expected that, if not detections, an 
interesting range of rate upper limits will be probed. 

A dedicated search for GWs from cosmic super-string cusps has been 
performed on S4 data. It has not revealed any GW signal and upper limits on 
the event rate have been set and used to exclude portions of the cosmic strings 
parameter space (string tension, loop size and reconnection probability). Big 
Bang Nucleosynthesis indirect exclusions are more stringent, but the GW 
search sensitivity obtained with more recent data is surpassing the indirect 
bounds in some portion of the parameter space [68] . 

3. Continuous wave searches 

3. 1 . Sources 

Neutron star are amazing objects for GW astronomy. Not only are they 
the most promising and well predictable source of GWs when they are part 
of a coalescing binary system, but spinning neutron stars can also radiate 
themselves a detectable amount of GWs. 

In the frequency band of ground-based GW detectors (10 Hz - 2 kHz), 
different emission mechanisms can lead to long-lasting, quasi-monochromatic 
GWs (see section HI of [22] and references therein). Among these, non 
axisymmetric deformations are considered the most reliable mechanism, at 
least for isolated systems. These deformations in the crust can be gener- 
ated by elastic stress [TOl [71] or by magnetic fields [72] . Mechanisms other 
than crustal distortions include GW driven instabilities of normal oscillation 
modes (r-mode ^73j and f-modes Chandrasekhar-Friedman-Schutz [71] insta- 
bilities), or free precession, when the star's rotational axis does not coincide 
with its symmetry axis. 

In the case of non-axisymmetric deformations, the gravitational wave fre- 
quency is twice the spin frequency and the amplitude of the emitted GW is 
proportional to the deformation which is expressed in terms of the elliptic- 
ity qj of the distorted object. Predictions for e are highly uncertain. The 
maximum deformation that a crust can support depends on the neutron star 
equation of state through the shear modulus of the inner crust and on the 
crust's structure (crystalline or amorphous), through its breaking strain [35] . 



6g _ ^Ex_ — w^ where the Is are the moments of inertia of the star and the spin axis is 
assumed to He in the z— direction. 
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Possible breaking strain values span several orders of magnitude, the largest 
ones (~ 10'^) corresponding to perfect crystals. Shear modulus values of the 
inner crust also span a few orders of magnitude, with solid strange-quark star 
models yielding the highest shear modulus values. Adding it all up maximum 
values of the ellipticity range from several 10~^ to a few 10~^. 

But what are the ellipticity values in actual pulsars ? This question does 
not have a clear-cut answer and it cannot be dismissed that a certain frac- 
tion of the neutron stars have negligible axial quadrupole moment. This is 
supported by observations of extremely low period derivates of millisecond 
pulsars that indicate that these recycled pulsars have little equatorial asym- 
metry (e < 10^^). Simulations of isolated spinning neutron star populations 
including gravitational and electromagnetic torques [75] provide insight in 
the type of constraints that even non-detections of continuous GW signals 
will place on neutron star population parameters in the advanced detector 
era. 

Accreting neutron stars in binary systems have a natural way to maintain 
large crust deformations, through the accretion spot. In fact it has been 
suspected for a long time that GW emission could balance the accretion 
torque and explain why spin frequency of known accreting neutron stars are 
never approaching the breaking limit [TOl [76]. However the actual search 
for signals from binary systems is significantly complicated by the additional 
unknown orbital parameters that one would need to search over [77] . 

3.2. The detection problem 

In the rest frame of the source and in absence of glitches or star quakes, 
continuous GWs have a constant amplitude and are quasi-monochromatic, 
with a slowly decreasing frequency (spin-down). Due to the relative motion 
between the source and the Earth, GWs are received at the detectors with a 
Doppler modulation. Additionally, due to the non uniform sensitivity of the 
detector across the sky, the signal amplitude also appears to be modulated. 
If one knows the frequency and frequency derivatives of the incoming GW 
signal and its origin in the sky, all the modulation effects can be removed 
and, broadly speaking, the detection problem broils down to that of detecting 
a sinusoidal signal, hence to integrating the data for long enough to accu- 
mulate sufficient signal-to-noise ratiqj This is conceptually what is done 



^The real situation is actually more complicated because of other unknown parameters 
(i.e. the inclination angle of the system, the polarization of the GW, the phase) but these 
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when targeting known objects for which the spin phase evolution, position 
and orbital parameters (if part of a binary) are known. 

For the vast majority of the known systems with GW frequency in the 
range of ground-based detectors we do not expect to be able to make detec- 
tions. This statement is based on comparing the typical strain (or amplitude) 
sensitivity of our search with the spin-down limit for the same system. The 
spin-down limit is given by the amplitude of a continuous GW that carries 
away energy at the observed rotational kinetic energy loss rate. It is the 
maximum GW amplitude that we can expect from a system for which we 
know the frequency evolution. For most known systems in the LIGO and 
Virgo frequency bands the spin-down limit is well below the values that our 
searches could detect. 

However there are many more neutron stars in our Galaxy than we know 
of (~ 10^), and the most interesting searches for continuous GWs aim at 
discovering a system that is spinning at the right frequency and is bumpy 
enough and close enough to be detected. The main challenge with this type 
of search is that one has to explicitly search over all parameters : position, 
frequency and spin-downs, if the system is isolated, and also over orbital pa- 
rameters for systems in binaries. The most sensitive techniques that combine 
the data coherently on longer and longer timescales T cannot be used because 
the resolution in parameter space increases rapidly with the observation time 
(oc T^ for searches including only first order spin-downs), thus making the 
problem computationally prohibitive. We hence resort to hierarchical meth- 
ods that iteratively identify at lower resolution regions in parameter space 
more likely to harbor a signal and then refine the search with higher reso- 
lution and sensitivity around such regions. The specific combination of the 
parameters that characterize the various stages is tuned depending on the 
data and on the parameter space that one wants to search. In the next section 
we discuss representative results of various types of searches for continuous 
wave signals. 

3.3. Results 

In the absence of a detection we place upper limits on the intrinsic GW 
strain amplitude Hq. Depending on the emission mechanism such upper limit 



may be maximized over or marginalized over and for the purpose of this discussion we can 
neglect them. 
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can be recast in terms of other physical quantities. If the emission is due to 
non-axisymmetric deformations then [78] 

levr^Gelf^ 

^0 = ^^— (1) 

where I is the moment of inertia, e the eUipticity, f the GW frequency and r 
the distance of the source. If one knows the distance to the putative source, 
as in the case of searches that target known objects, it is straightforward to 
see how upper limits on ho may be recast as upper limits on the ellipticity 
at a fixed nominal value of I (typically chosen 10^*^ kgm^ ) or as an exclusion 
region in the I — e plane, allowing I to vary. 

When the upper limits stem from a search that does not target any specific 
object or an object whose distance is unknown, the ho upper limits can be 
recast as lower limits on the distance of the source having assumed that it 
is emitting GWs because of a non-axisymmetric deformation and that it is 
doing so at the spin-down limit for some value of the spin-down. 

Under the assumption that the GW emission stems from r-modes, Eqs 
(23) and (24) of [79j allow to convert ho upper limits in upper limits on the 
dimensionless r-mode amplitude a for fiducial values of the source parameters 
such as mass, radius and equation of statq^ 

Targeted searches (known position and frequency): The ephemerides of 
many pulsars are regularly updated thanks to optical. X-ray and 7- 
ray pulsar observations. Such ephemeris are then used to determine 
the putative signals' Doppler and amplitude modulations employed by 
all methods developed to search for known pulsars in GW data: i) 
a matched-filter method [75] [HO] ii) a time-domain heterodyne search 
which generates Bayesian posteriors on the parameters associated with 
the signal (amplitude, source orientation, initial phase and polariza- 
tion angle) [81] iii) a Fourier-domain determination of the signal at 
the 5 frequencies (carrier plus 2 pairs of sidebands) at which the sig- 
nal power is split by the amplitude and phase modulation due to the 
sidereal variation of the detector's beam pattern functions [HS]. The 
first continuous wave search performed with LIGO data targeted one 



^Note however that the conversion from ho upper hmits derived from detectabihty 
studies of a specific search on a specific population of signals is correct if the polarization 
angles of the injected population were uniformly distributed. 

17 



pulsar, J1939+2134 {frot = 64:2Hz) using data of 2002 [H3]. Since then, 
the number of known pulsars searched for in each new data set has in- 
creased steadily while more pipelines have been developed and refined. 
This, combined with more sensitive data and longer data sets has led to 
more constraining upper limits over the time [HHIHSllHn]- With S5 data 
the lowest strain amplitude upper limit is 2.3 x 10^^^ achieved for PSR 
J1603-7202 and the lowest limit on the ellipticity is 7.0 x 10"^ for PSR 
J2124— 3358. The ho upper limits values from our GW observations are 
more constraining than the spin-down limits for two young pulsars: the 
Crab pulsar (/rot=29.78 Hz) and the Vela pulsar (/roi=11.2 Hz). For 
the Crab the 95% confidence upper limit on the strain amplitude is 
2.0 X 10~^^, yielding a 2% upper limit value on the fractional rota- 
tional energy loss through GW emission [871 [86] . For Vela the strain 
amplitude limit is 2 x 10~^^, bounding to less than 35% the amount 
of rotational energy j88] lost in GWs. Albeit the non-detection per 
se is not unexpected since the ellipticity values corresponding to the 
spin down limits are rather high (a few 10~^ and 10~^ for the Crab and 
Vela pulsars respectively), the constraints on the fraction of energy loss 
in GWs are informative measurements that could not be obtained by 
other means. 

Directed searches (known position but unknown frequency): These searches 
are somewhere in-between the targeted searches described above and 
the all-sky searches for completely unknown objects that will be de- 
scribed next. By directed searches we mean searches aimed at a given 
direction in the sky where a compact object, most of the time a non 
pulsing neutron star, may be emitting strong continuous GWs. This 
is for instance the case of the supernova remnant Cassiopeia A (~ 3.4 
kpc) which is suspected to host a compact central object born less than 
400 years ago. Under the assumption that its spin frequency evolution 
has been largely governed by GW emission, an indirect limit on the 
GW strain amplitude, h.^^^, can be derived which only depends on the 
age of the object 



h;«. . 1.2 X io-(^^)^(^^)(?52Z5^) (2) 

with r being the age of the object [89j. In the absence of any ob- 
servational information on the rotation period of the central object, a 



200 Hz wide frequency band was searched for continuous signals from 
this object in S5 data [90]. The search band was chosen so that the 
minimum detectable strain amplitude across that band would be lower 
than the hl^^^ limit. No significant GW candidate was identified and 
95% confidence upper limits on the strain amplitude (7 x 10~^^ at 150 
Hz), ellipticity (4 x 10"^ at 100 Hz), and r-mode amplitude (5 x 10"^ 
at 300 Hz) were derived l90l. 



The low mass X-ray binary Sco X-1 (2.8 kpc), the brightest X-ray 
source in the Galaxy, is expected to be also a strong GW emitter. The 
detection of continuous GWs from this system would allow to test the 
torque balance mechanism of accreting neutron stars. The rotation 
frequency of Sco X-1 is unknown making the search over frequency, 
spin-down and orbital parameters computationally very challenging. A 
20 Hz wide frequency band was searched in LIGO 2003 data using of 
only few hours of the most sensitive data. The upper limit on the strain 
amplitude is few orders of magnitude higher than the torque balance 
prediction, which is ~ 3 x 10~^^ for a source emitting GW at 540 Hz 



All-sky, all-frequency searches (objects not previously known): 

As described above, all-sky, all-frequency searches for unknown neutron 
stars are the most challenging of the continuous GW searches carried 
out to date. These searches typically do not begin before all the data to 
be analyzed has been acquired. The core search lasts of the order of a 
year and the post processing takes a comparable amount of time. This 
means that results from this type of analyses lag behind the results for 
other types of signals. The most recent results come from two searches 
on S5 data [211 E2]- In [SI] we add incoherently 30-minutes long power 
spectra over the entire duration of the run, spanning a frequency range 
between 50 and 800 Hz. This is a fast-turn-around of results type 
of search. In ^2] the searched frequency range extends to 1200 Hz 
(the higher frequencies are the most computationally expensive), one 
year of S5 data is considered, split up in 121 30-hour segments where 
the data from the two detectors is coherently combined. This is a 
very compute-intensive search and it leveraged the compute power of 



19 



Einstein@Honiqj, a volunteer computing project where members of the 
general public donate compute cycles from their personal computers in 
order to help with a scientific endeavor. Both methods (Powerflux and 
Einstein@Home) set upper limits on the strain amplitude as function of 
the GW frequency which are comparable. For instance, the best 90% 
CL upper limit achieved by Einstein@Home is 7.6 x 10~^^ at 152.5 Hz 
in agreement with the expected sensitivity of the analysis. The reach 
of the search extends to a few kpc around 150 Hz. 

The most recent results of continuous wave searches are summarized in 
Figure |3| 

4. Stochastic gravitational wave background searches 

4.I. Sources 

A cosmological and astrophysical stochastic background of GWs (SGWB) 
is expected from the incoherent superposition of GWs emitted by many 
sources. Cosmological sources include amplification of quantum vacuum fluc- 
tuations during inflation [93j, electro- weak phase transitions [91], pre Big 
Bang models [95], [96| and cosmic (super-) strings [971 EHl |99]. Most of these 
processes took place in the early stages of the universe making GWs a unique 
way to probe primordial evolution. An astrophysical SGWB may be gener- 
ated by the concurrent emission of a large population of sources such like 
compact binary coalescence and merger, core collapse supernovae or rotating 
neutron stars. 

A SGWB is described in terms of the GW energy density 

ilawif) = i^ (3) 

Pc d\nf 

where pc is the critical energy density of the universe and dpcw/dlrif is the 
GW energy density per logarithmic frequency interval. Predictions on flow 
as well as indirect phenomenological and experimental bounds are illustrated 
in figure 111 which is taken from |100] . The most stringent bound comes from 
Big Bang nucleosynthesis (BBN) and cosmic microwave background (CMB) 



'http : //einstein . phys . uwm . edu . 
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Figure 3: The indirect spin-down upper limits (open circles) for the known pulsars are 
below the measured values (open triangles) for all pulsars but the Crab (at about 60 Hz) 
and the Vela pulsar (at about 22 Hz). The ft-o upper limits from these searches are the 
lowest compared to the ones from other searches because the parameter space is much 
smaller : virtually one template versus 10^^ or 10^* templates for the Cassiopeia A (black 
small stars) and the EinsteinQHome (dark grey crosses) searches, respectively. This affects 
the sensitivity in two ways 1) through the trial factor, i.e. the largest random fluctuations 
increase over a larger set of independent trials and this is why a signal needs to be larger 
in order to be detected with the same confidence over a larger parameter space 2) the most 
sensitive and computationally intensive coherent search techniques can be employed when 
the same search has to be repeated fewer times. The Powerflux (light grey open squares) 
search examined more than 10* points in parameter space but with a non coherent method. 
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measurements, and it depends on number N^ of neutrino families at the time 
of the BBN [M] : 



^BBN = / nGw{f)d\nf < 1.1 X 10-^(iV, - 3) 



(4) 



Primordial nucleosynthesis provides an upper bound A'^^ — 3 < 1.4. The CMB 
data can be used to derive a limit on the maximal SGWB energy density at 
the time of photon decoupling |102j 



^CMB = / nGw{f)dlnf < 1.3 X 10' 



(5) 



These are the two phenomenological upper bounds that can be directly com- 
pared to direct LIGO- Virgo measurements. 
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Figure 4: Comparisons of different SGWB search upper limits on tlie GW spectrum ^gw 
and indirect bounds IIOOI. 
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4-2. Results 

Since the signal is a stochastic process the search is made by cross- 
correlating detector strain amplitudes (the GW amplitude is much smaller 
than the instrumental noise and the signal is expected to be continuous). For 
cosmological SGWB sources one assumes that the SGWB is isotropic and en- 
tirely described by a spectral power density whose frequency dependence is 
parameterized by a power law with index a, Qcwif) = ^aif/fr)"", where fr 
is an arbitrary reference frequency. On the contrary for astrophysical sources 
it makes sense to search a background with a spatial distribution, as if it was 
coming from a point in the sky. 

The sensitivity of each pair of detectors depends on the noise power spec- 
trum of each detector and on the observing geometry (detectors location and 
orientation and spatial distribution of the background). Having co-located 
and co-aligned detectors like the two Hanford LIGO detectors maximizes the 
search sensitivity over the full frequency band of the detectors, but taking 
into account correctly the presence of correlated noise that could cancel, cor- 
rupt or mime the measurement of a SGWB is a serious issue. On the contrary 
the LIGO Hanford and Livingston pair presumably do not have correlated 
noise and benefit from their good relative alignment that yields good sen- 
sitivity in the low frequency region where the detectors' sensitivity is the 
highest. The isotropic search carried out in the 4L5 — 169.25 Hz frequency 
band with S5 data of this pair of detectors set a direct upper limit on Qq at 
95% confidence level of 6.9 x lO'^ [IDQ]. A SGWB of this strength would 
contribute 9.7 x 10~^ to the total energy density. This limit hence beats 
the BBN and CMB bounds. It also constrains values of the cosmic strings 
free parameters (string tension fi and loop size e) not excluded by any other 
experiments; for a rather low value of the string reconnection probability 
10"^, 7 X 10"^ < Gfx < 1.5 X 10"'^ for e < 8 X 10"^^ Measurements from 
Virgo, LIGO Hanford and LIGO Livingston probe the existence of a SGWB 
at higher frequency where the Hanford-Livingston pair sensitivity is slightly 
worse. An upper limit on ^^3 < 0.33 at 900 Hz has been set combining LIGO 
Hanford, LIGO and Virgo pairs |103j . This improves a limit set by LIGO 
Livingston and the ALLEGRO bar detector (fiGw(915Hz) < 1.02) [TM] . 

Finally, searches for point-like or more extended sources have been con- 
ducted over the same data, producing upper limits on the strain amplitude in 
each 0.25 Hz wide frequency bin. Three specific directions (the nearby low- 
mass X-ray binary Sco X-1, the galactic center and SN 1987A) have been 
studied, but no indication of a SGWB has been found |105j . 
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5. Towards the second generation of detectors 

In a few years from now a second generation of GW detectors will come 
online with initial improvements in sensitivity of about a factor of two with 
respect the the S6/VSR2-3 runs and ramping up over the space of a few years 
to sensitivities about of factor of ten higher fS]. The network will include 
two LIGO detectors at Hanford and Livingston, USA, a Virgo detector in 
Cascina, Italy, and the smaller scale GEO detector in Germany. By the 
beginning of the next decade one of the LIGO detectors is likely to become 
operative in India and soon thereafter a large scale detector is expected in 
Japan enhancing the sky coverage and source localization capabilities of the 
global network. 

The science runs of the second generation detectors promise to yield the 
first direct detection of GWs by 2020. With an order of magnitude better 
sensitivity, and despite neutron star merger rate uncertainties, it is likely 
that the first astrophysical events observed be one of the 1 - 400 neutron 
star mergers predicted per year [Hj. Search pipelines aimed at detecting 
such signals have been exercised in blind injection challenges [51] and have 
demonstrated their readiness to confidently detect such events. However 
there may be surprises and the broad variety of GW searches that have 
been developed in the past years reflects the breadth of the net that we 
cast. With the flrst rapid EM follow-ups of GW triggers and with the GW 
searches triggered by EM and neutrino observations we have taken the flrst 
steps towards multi-messenger astronomy which will be the new paradigm 
for the study of compact objects soon after the flrst few GW detections. A 
very exciting decade lies ahead ! 
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